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Abstract. Higii spatial resolution observations in the 1 to 3.5 fim region of the Galactic Center source known 
historically as IRS 13 are presented. They include ground-based adaptive optics images in the H, Kp (2.12/0.4 ^m) 
and L bands, HST-NICMOS data in filters between 1.1 and 2.2 /xm, and integral field spectroscopic data from 
BEAR, an Imaging FTS, in the He I 2.06 /xm and the Br7 line regions. Analysis of all these data provides a 
completely new picture of the main component, IRS 13E, which appears as a cluster of seven individual stars 
within a projected diameter of ~0.5" (0.02 pc). The brightest sources, 13E1, 13E2, 13E3 which is detected as 
a binary, and 13E4, are all massive stars of difi^erent type. The star 13E1 is a luminous, blue object, with no 
detected emission line. 13E2 and 13E4 are two hot, high-mass emission line stars, 13E2 being at the WR stage 
and 13E4 a massive O-type star. In contrast, 13E3A and B are extremely red objects, proposed as other examples 
of dusty WR stars, like IRS 21 (Tanner et al. 2002). All these sources have a common westward proper motion 
(Ott et al. 2003) indicating they are bounded. Two other sources, detected after deconvolution of the AO images 
in the H and Kp bands, are also identified. One, that we call 13E5, is a red source similar to 13E3A and B, 
while the other one, 13E6, is probably a main sequence O star in front of the cluster. Considering this exceptional 
concentration of comoving massive hot stars, IRS 13E is proposed as the remaining core of a massive star cluster, 
which could harbor an intermediate-mass black hole (IMBH) (Portegies Zwart & McMillan 2002) of ~13OOM0. 
This detection plays in favor of a scenario, first suggested by Gerhard (2001), in which the helium stars and the 
other hot stars in the central parsec originate from the stripping of a massive cluster formed several tens of pc 
from the center. This cluster would have spiraled towards SgrA*, and IRS 13E would be its remnant. Furthermore, 
IRS 13E might be the second black hole needed according to a model by Hansen & Milosavljevic (2003) to drag 
massive main-sequence stars, in the required timescale, very close to the massive black hole. The detection of a 
discrete X-ray emission (Baganoff et al. 2003) at the IRS 13 position (within the positional accuracy) is examined 
in this context. 

Key words, instrumentation: adaptive optics - instrumentation: NICMOS - infrared: stars - X-ray: source - 
Galaxy: Center - stars: Wolf-Rayet 



1. Introduction 

In the early images of the central region of the GC 
recorded in the near infrared, at the best seeing-limitcd 
resolution, several bright point sources dominate the 
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--20"x20" field centered on SgrA*. With the radical im- 
provement of angular resolution through multiple short 
exposures with shift-and-add (SHARP camera, Eckart et 
al. 1995) or speckle techniques (Ghez et al. 1999), with 
the images from the NICMOS cameras on board HST 
(Stolovy et al. 1999), to the advent of adaptive optics 
correctors behind large telescopes as NAOS (Ott et al. 
2003), a more complex vision of this crowded field has 
emerged. At a spatial resolution in the best case of 0.05" 
(Ghez et al. 1999) new, faint stellar sources appear and the 
early-identified sources are often resolved into several com- 
ponents. Therefore, the observed spectra, generally made 
at a resolution limited at best by a slit not less than 1" 
wide, are actually composites of emission from stellar ob- 
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jects of different spectral type as well as from emission 
lines from the surrounding interstellar medium (ISM). The 
stellar type identified from these spectra can be wrong if 
attributed to a single source, which can erroneously look 
unusual. This care led to the paper on the revised iden- 
tifications of the central cluster of He i stars by Paumard 
et al. (2001) made from slitless integral field spectroscopy, 
with BEAR an imaging Fourier Transform Spectrometer 
(Maillard 2000). This instrument makes it possible to ob- 
tain near-infrared spectroscopy at the seeing-limited reso- 
lution of Mauna Kea (i.e. ~0.5"), which represents a sig- 
nificant improvement. However, this resolution is not suf- 
ficient for such a crowded field. Therefore, the data were 
compared to adaptive optics (AO) images in the K band 
of the same field, with a spatial resolution of ^^0.15", in 
order to check whether the emission sources at the angular 
resolution of seeing were single stars or not. 

Among the sources studied in this paper the object 
historically named IRS 13, is a typical example. Located 
approximately 3.6" south-west of SgrA*, at the edge of 
the circular structure of the SgrA West HIT region called 
the Minicavity, it appears in early works at near-infrared 
wavelengths, in J, H, K (Rieke et al. 1989) and L (Allen 
& Sanders 1986), as a bright spot. However, from a dis- 
cussion on the sources of energy at the Galactic Center, 
Rieke et al. (1989) note "there are separate luminosity 
sources in the core of the 10 sources 1, 9 and 13", 
because, particularly for IRS 13, the energy distribution 
from 1 to 5 /xm presents a sudden, steep increase beyond 
3 /im. A first work at subarcsecond resolution, by lunar 
occultation in the K band (Simon et al. 1990), indicates 
that IRS 13 resolves into a pair of sources, separated by 
~ 1.2", which were henceforth designated as IRS 13E and 
13W. In the photometric survey of Ott ct al. (1999) from 
SHARP imaging data (Eckart ct al. 1995), at a limit of res- 
olution of ^^0.15" obtained by deconvolution, two equally 
bright components (Kmag = f 0.26) arc reported, IRS 13E1 
and 13E2, with a separation of 0.2". Paumard ct al. 2001) 
published the first AO map of IRS 13E in the K band, ex- 
tracted from a larger image of the central region obtained 
with the CFHT-AO system (Lai ct al. 1997). We noted 
that a fainter third source, we called IRS 13E3, was also 
present, forming a kind of equilateral triangle with 13E1 
and 13E2. In parallel, Clenet et al. (2001) published a pho- 
tometric analysis of the same data. For IRS 13E they re- 
ported three components in the K band, noted 13E1, 13E2 
and 13N. From the given offset coordinates the source 
called 13N in this paper is not exactly coincident with 
the source we had previously called 13E3 (Paumard et al. 
2001). 

In the meantime, spectroscopic works on the stellar 
population of the central parsec were conducted. Several 
spectra of IRS 13 (Blum et al. 1995, Libonatc ct al. 1995, 
Tamblyn et al. 1996), and specifically of IRS 13E ( Krabbe 
et al. 1995, Genzel et al. 1996, Najarro et al. 1997), and 
of IRS 13W (Krabbe ct al. 1995) have been published. 
They cover mostly the K band, i.e. all or part of the 1.95 
- 2.45 /um range, and one the 1.57 - 1.75 /xm region of 



the H band (Libonate et al. 1995). From these spectra 
IRS 13W is unambiguously a cool star with the strong vib- 
rot CO signatures at 2.3 /im. Emission lines were detected 
on IRS 13E, typical of a luminous, hot star: strong Hei 
2.058, 2.112 /im, Br7 line and other Brackett lines up to 
Brl2, plus weaker lines of [Fell], [Fciii], and a weak emis- 
sion at 2.189 /im attributed to Heii. The 7-10 line is the 
most intense of this ion in the K band, and is supposed un- 
blended, while the 8-14 line cannot be distinguished from 
Br7 of which it is separated by less than 2 cm^^. The He ii 
line if detected is a precious type indicator in the classifica- 
tion of WR stars (Tamblyn et al. 1996, Figer et al. 1997). 
From these spectral characteristics Libonate et al. (1995) 
concluded "the IRS 13 spectrum bears a strong resem- 
blance to the low-resolution K-band spectra of P Cygni 
(a LBV) and the AF source". Krabbe et al. (1995) pro- 
posed IRS 13E as a WN9 star. From a K-band spectrum of 
the Minicavity Lutz et al. (1993) attributed to [Feiii] sev- 
eral emission lines (2.145, 2.218, 2.242, 2.348 /xm). Images 
at 2.218 /xm show that this line is present all over the 
Minicavity, with a particular enhancement at the posi- 
tion of IRS 13. From NICMOS data in the F164N filter 
(Stolovy et al. 1999) the lower ionized Fe species, [Fell] 
is detected in emission in the central ISM, particularly 
strong at the edge of the Minicavity, but remarkably ab- 
sent at the position of IRS 13E. As mentioned in this 
paper the ionization condition for [Feiii] requires 16.2 eV 
while only 7.2 eV are required for [Fell]. Hard ionization 
radiation is likely originating from IRS 13E. Hence, it can 
concluded that these iron lines are not of stellar origin but 
are present in the surrounding ionized gas. 

From the extraction of the He i 2.058 /xm line profile at 
high spectral resolution from the BEAR data, Paumard et 
al. (2001) concluded that one of the three sources identi- 
fied as forming IRS 13E, instead of a LBV- type star should 
be already at the WR stage. The main argument was 
the width of the observed emission line profile (FWHM 
~ 974km s~^), typical of WR stars, making the source 
belonging to a class of broad-lino stars including 8 other 
stars in the central cluster of helium stars. It was also mea- 
sured that the broad-line stars were in average weaker by 
~2.4 mag in the K band than the rest of helium sources 
which were characterized by a narrow Hei emission line. 
With a Kmag for IRS 13E3 consistent with the magni- 
tude measured for the other broad-line stars, and 13E1 
and 13E2 much brighter, it was proposed that IRS 13E3 
should be the WR-type helium emitter, of because and 
It was not possible to precise the 13E1 and 13E2 stellar 
type. 

In the centimetric domain, Zhao & Goss (1998) pre- 
sented the detection of IRS 13 at 7 and 13 mm with the 
VLA, at a resolution of 0.06", which they reported to be 
the brightest radio continuum source after SgrA* at the 
Galactic Center. They resolved the source into two com- 
ponents, one with no significant proper motion while the 
other one is moving south at a rate 6.2 ±1.1 mas yr~^. 
They called the two compact Hii regions IRS 13E and 
IRS 13W, which was improper, since these denomina- 
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tions had already been given to infrared sources at IRS 13 
as reminded above, with which they are not coincident. 
However, this detection is another element which makes 
this source special. 

The detection of a bright, discrete X-ray emission 
source (Baganoff ct al. 2003, Muno ct al. 2003) at the 
IRS 13 position within the positional accuracy - source 
CXOGC J174539.7-290029- among the brightest sources 
within the central parscc besides SgrA*, is a last clement 
contributing to make this source an object of interest. 
This coincidence already triggered the interest of Coker 
et al. (2002) who presented the first Chandra X-ray spec- 
trum of IRS 13 and deduced that it was consistent with 
a highly absorbed X-ray binary system. They concluded 
that IRS 13E2 should be a compact post-LBV binary 
whose colliding winds were the source of the X-ray emis- 
sion. 

As a conclusion, the origin of the brightness of IRS 13, 
one of the brightest objects at all wavelengths in the vicin- 
ity of SgrA*, has been a matter of debate. We examined 
all the high-angular resolution images in the near infrared 
currently available. This analysis made it possible to build 
a completely new picture of IRS 13 which is presented in 
this paper. We will show that the peculiar spectral energy 
distribution (SED) previously reported is well explained 
by the nature of the individual sources which compose 
IRS 13. Unexpectedly, this analysis provides new insight 
to the unsolved question of the formation of the central 
massive cluster. It has also given the opportunity through 
the study of the stars in the IRS 13 field to character- 
ize a sample of the star population in the central parsec. 
A preliminary version of the paper was presented at the 
Galactic Center Workshop 2002 (Maillard et al. 2003). 

2. High-angular resolution images of IRS 13 and 
data reduction 

Ground-based AO images from several telescopes and 
space-based NICMOS data containing IRS 13 in their 
field have been gathered. Their characteristics are given 
in Table 1. 

2.1. Adaptive Optics data 

Wc have analyzed data from three different AO systems: 
PUEO/CFHT (Lai ct al. 1997), Hokupa'a/Gemini North 
(Graves et al. 1998), and Adonis/ESO 3.6-m telescope 
(Beuzit et al. 1997). The CFHT data, in the K band, were 
obtained on 26 June, 1998. They were already presented 
and analyzed in detail in Paumard et al. (2001). The to- 
tal field covers 40"x40" centered on SgrA*. The FWHM 
of the Point Spread Function (PSF) varies from 0.13" to 
0.20" in the field since PUEO has a visible wavefront sen- 
sor, requiring to use a V ~ 14 star 20" to the north-east 
of the field center. The L-band data were obtained with 
the ADONIS visible wavefront sensor in 2000, from 20 to 
22 May, on the ESO 3.6-m telescope. These data are de- 
scribed in Clenet et al. (2001). The final L-band image 



we used has a PSF with a FWHM of 0.291" and covers a 
field of ~ 13"xl3" centered on SgrA*. The Gemini North 
data were part of the AO demonstration nm conducted 
by F. Rigaut on the Galactic Center in July 2000. The 
data were obtained with the Hokupa'a AO system and 
the QUIRC camera (Graves et al. 1998) in the Kp and 
the H band (Table 1) respectively 3 and 6 July in field 
1, centered on SgrA*. The field coverage of each image is 
20"x20". For the H image the FWHM of the PSF varies 
from 0.115 to 0.19" and for the Kp image from 0.12 to 
0.18". In the vicinity of IRS 13 the measured FWHMs are 
respectively 0.180" and 0.172". These AO data are char- 
acterized by a low Strehl ratio of 2.5% in H and 7% in 
Kp. The small portion of the AO image in the Kp filter 
analyzed in the paper is shown on Fig. 1. 

2.2. NICMOS data 

Six filters, coded FllOM, F145M, F160W, F187N, F190N, 
and F222M were used in observing the stars in the in- 
ner parsec of the Galactic Center with the NICMOS cam- 
eras on board HST, as part of three independent Galactic 
Center programs (7214, 7225, and 7842). These data were 
taken between Aug. 1997 and Aug. 1998. The filter proper- 
ties (central wavelength, bandwidth, zero-magnitude flux, 
pixel size) and the data processing of the raw data arc de- 
scribed in Rieke (1999) and Stolovy et al. (1999). Although 
each program covered a larger region, we extracted a small 
portion of ~ 4" x 4" centered on IRS 13E of the image from 
each filter. All these diffraction-limited images were partic- 
ularly useful to derive the SED of the stellar components 
of IRS 13 and its environment. 

Apart from wide (W) or medium (M) bandpass fil- 
ters, F187N is a 1% narrow-band filter centered on the 
Paa line. By subtraction of F190N, a comparable narrow- 
band filter in the nearby continuum, the distribution of 
the ionized gas in the inner parsecs was obtained (Stolovy 
et al. 1999). In this map are also discrete emission sources 
which must come from Paa emission in the atmosphere of 
hot stars. We used a best-fit scale factor between F187N 
and F190N in order to minimize both positive and neg- 
ative stellar residuals of the stars in the central parsec. 
Negative residuals can arise in stars with Paa absorption 
and/or in stars with local excess extinction along the line 
of sight. For this paper, we identified emission line stars as 
those with very significant F187N/F190N ratios (exceed- 
ing 1.2). The portion of the resulting Paa image with the 
IRS 13 complex was essential in determining which of the 
individual stars was an emission line star (Fig. 3). 

2.3. Data reduction 

The initial processing of each dataset is not described here. 
It can be found in the papers cited with the presentation of 
each of them. Below are listed the main operations which 
were performed subsequently to the initial data reduction. 
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Table 1. Instruments, filter properties and year of data acqui- 
sition of the high-resolution images of the IRS 13 field 



Filter Instrument / Telescope A/AA year 

(/jm) 

FllOM NICMOS/HST 1.10/0.200 97-98 

F145M NICMOS/HST 1.45/0.197 97-98 

F160W NICMOS/HST 1.60/0.400 97-98 

H Hokupa'a-hQuirc/Gcmini N 1.65/0.296 2000 

Paa NICMOS/HST 1.87/0.019 97-98 

F190N NICMOS/HST 1.90/0.017 97-98 

Kp Hokupa'a-hQuirc/Gemini N 2.12/0.410 2000 

K PUEO-hKIR/CFHT 2.20/0.340 1998 

F222M NICMOS/HST 2.22/0.143 97-98 

L ADONIS-|-COMIC/3.6m ESO 3.48/0.590 2000 



not applied to the CFHT K-band nor to the NICMOS im- 
ages, which both PSFs show significant secondary rings, 
because they arc not correctly handled by this code, writ- 
ten for seeing-limited images. Anyway, it was applied on 
the Gemini/Hokupa'a data because of their low Strehl ra- 
tio. 

2.3.2. Calibration 

All the NICMOS data were calibrated as described by 
Rieke (1999). The calibration of the sources in the L band 
was based on the photometry of IRS 13W reported by 
Clcnet et al. (2001), a prominent, isolated source in this 
band. The Gemini data in the H and the Kp filters were 
not calibrated. In order to minimize photometric uncer- 
tainties between data from various origins we calibrated 
the two filters by interpolation from the NICMOS data 
which offer a set of filters close to H and Kp. Based on 
IRS 13E1, a bright and wcll-mcasurcd star by MCS dccon- 
volution in the stellar complex, and taking into account 
the exact filter bandpasses, we applied a strict linear in- 
terpolation between the F160W and F190N flux of this 
star for the H band (central wavelength 1.65 /xm), between 
F190N and F222M for the Kp band (central wavelength 
2.12 /im, sec Table 1). With this source calibrated that 
way in the H and the Kp filters the flux measurements 
of all the other sources in the two filters were calibrated. 
The intensities of the detected sources in the 8 filters were 
expressed in /xJy as were originally the NICMOS data. 

2.3.3. Proper motions 

We obtained from Ott et al. (2003), who have conducted 
an analysis of ten years of SHARP data, providing more 
than 1000 proper motions in the central parsec, the proper 
motions of the main components of IRS 13 and of its 
nearby field. 

2.4. BEAR data 

Even though the BEAR data are not at the spatial reso- 
lution of the NICMOS and AO data, the Br7 and 2.06 /xm 
He I line profile at IRS 13E were used as complementary 
information to help determine the spectral type of the un- 
derlying stars. The IRS 13 complex is located in a region 
of interstellar emission, intense in Br7 (Morris & Maillard 
2000, Paumard et al. 2003a) as well as in the 2.06 /xm he- 
lium line (Paumard et al. 2001). Hence, the line profiles 
had to be corrected to remove the ISM contribution, made 
here of two main velocity components (Fig. 4 - central 
panel), superimposed to the stellar profile. This operation 
required special attention for the Br7 profile since the ISM 
emission is locally very intense. Only the high resolution 
(21.3 km s"""^) made it possible to separate, with some ap- 
proximation anyway, the stellar profile which is naturally 
broader, from the interstellar components. 



2.3.1. Methods of star extraction 

To perform the extraction of the individual stars from the 
AO and NICMOS data, StarFinder, an IDL procedure 
(Diolaiti et al. 2000) specially written for AO data, was 
first used. The preliminary operation consists of building a 
good PSF by averaging several isolated and bright stars in 
each image. By adjusting the PSF to the local peaks in the 
field the exact position and the flux of the corresponding 
stars are retrieved by the procedure. 

The image residuals left by applying this procedure led 
us to suspect that more sources might be present in the 
IRS 13 complex, which could not be detected because the 
spatial resolution was not high enough to separate them. 
To improve the star detection we applied a new decon- 
volution code called MCS (Magain et al. 1998). Contrary 
to the StarFinder procedure, the MCS program uses an 
analytical PSF (a Moffat- type function). The final PSF 
is chosen prior to deconvolution, with a width compatible 
with the original image sampling. Thus, the final PSF can 
be narrower than the observed one without violating the 
Shannon sampling theorcuri. The contribution of a contin- 
uous background is also niatclicd to the image with an 
adjustable smoothing. The sampling was high enough for 
the Gemini AO images, in the H and the Kp bands, to 
provide a substantial gain in resolution. Because of the 
sampling of the L-band image a more limited gain was 
obtained, which was useful anyway. For the H and the Kp 
images the width of the synthetic PSF is equal to 0.040" 
and to 0.192" for the L-band image, i.e. a gain in resolu- 
tion respectively of a factor 4.5 in H, 4.3 in Kp and 1.5 
in L. However, this method applied to the IRS 13 images 
demanded many trials to reach a stable and plausible solu- 
tion because of the presence of a non uniform background 
in which the sources are embedded. The convergence was 
helped by the comparison between the solutions for the 
H and Kp filters. We imposed the sources to be detected 
in both bands, to avoid the risk of taking deconvolution 
artifacts for sources. The two filter bandpasses are almost 
adjacent which makes pertinent the detection in both fil- 
ters, even for very red or very blue objects. This code was 
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3. Results 

From the 8 high-resolution images between 1.1 and 3.5 /xm 
(Table 1) and the spectroscopic data available, the follow- 
ing results on the IRS 13 complex and its environment 
have been obtained: 

3.1. Star detection 

Twenty sources in total (Table 2) are detected both in 
H and Kp after MCS dcconvolution of the Gemini data 
in the small IRS 13 field, which are identified on Fig. 2. 
All these sources were searched for with StarFinder for 
the filters where the dcconvolution operation could not 
be applied, and again with the MCS code in the L band. 
Empty positions in Table 2 indicate that the source is 
not detectable in a specific filter. An upper limit is given 
at some positions for the extreme filters (1.1 /xm and L) 
where the detection is the most difficult. This limit is not 
uniform in the field depending on the proximity of another 
bright source, particularly in L. 

Table 2. Photometric measurements in all the filters, in 

log(/iJy), of all the stars detected by MCS deconvolution 
(Sect. 2.3) of the 2.5" x 2.5" IRS 13 field from the H and Kp 
Gemini AO images. Measurements in the NICMOS filters (F) 
were obtained by StarFinder, and in the L band by MCS anal- 
ysis. The sources are listed by order of decreasing brightness 
in the Kp band, except IRS 13W, which is listed first. Values 
preceded by < in the 1.1 /xm and the L bands are minima of 
detection, i.e. flux upper limits. 



ID 


FllOM F145M F160W H F190N 


Kp 


F222M 


L 


W 


<0.56 


2.58 


3.28 


3.20 


3.96 


4.27 


4.36 


4.89 


El 


1.76 


3.40 


3.82 


3.93 


4.27 


4.40 


4.49 


5.03 


E2 


1.61 


3.22 


3.71 


3.80 


4.17 


4.38 


4.48 


5.37 


E4 


0.93 


2.68 


3.17 


3.28 


3.69 


4.11 


4.18 


<4.66 


5 


0.78 


2.66 


3.20 


3.30 


3.72 


4.04 


3.93 




6 


0.39 


2.51 


3.07 


3.14 


3.57 


3.94 


3.89 




7 


1.23 


2.78 


3.20 


3.26 


3.63 


3.92 


3.88 




8 


0.62 


2.55 


3.08 


3.08 


3.53 


3.90 


3.81 




E3A 


<0.45 




2.29 


2.43 


3.11 


3.81 




5.46 


10 


0.11 


1.88 


2.63 


2.81 


3.23 


3.68 


3.55 




E3B 


<0.45 




2.11 


2.24 




3.57 




5.29 


12 


0.69 


1.98 


2.62 


2.66 


3.15 


3.48 


3.45 




13 


-0.09 


2.02 


2.58 


2.65 


3.14 


3.46 


3.40 




14 


-0.16 


2.01 


2.55 


2.65 


3.11 


3.43 


3.34 




15 


-0.06 


1.75 


2.25 


2.37 


2.82 


3.12 


3.23 




16 




1.57 


2.18 


2.36 


2.71 


3.11 


2.92 




E5 


<-0.17 


1.60 




1.75 


2.72 


3.08 


3.61 


5.07 


18 


0.56 


1.87 


2.27 


2.30 


2.74 


2.99 


3.02 




E6 


<0.48 




2.46 


2.53 




2.97 




<4.68 


20 


0.86 


1.42 


2.01 


2.11 


2.57 


2.88 


2.89 





In conclusion, IRS 13E is resolved into a compact clus- 
ter of at least 7 objects encircled within ~0.5" (Fig. 2). 
The two brightest sources, 13E1 and 13E2, had already 
been identified in Ott et al. (1999) and in Paumard et al. 



(2001). The one we had noted 13E3 appears double after 
deconvolution. Thus, we call the two components 13E3A 
and 13E3B. By continuity, the closest bright source, north 
of 13E3, is called 13E4. Two other sources revealed by de- 
convolution are called 13E5 and 13E6. 

3.2. Astrometry 

The positions of the sources in the 2.5" x 2.5" IRS 13 field 
in offsets from SgrA* are given in Table 3. The astrometry 
was retrieved from the results of the deconvolution of the 
portion of the Gemini Kp-band image. The precision of the 
positions depends on the brightness of the sources. For the 
brightest sources the relative 1-a position error is equal to 
1.0 mas. From these positions the projected angular sepa- 
rations between the main IRS 13E components are given 
in Table 3.2, translated to AU, by taking a distance to the 
GC of 8 kpc (Reid 1993). All the sources of the small field 
are detected by the reprocessing of the SHARP data we 
obtained prior to publication from Ott et al. (2003), except 
that a single source is given for E3A and E3B (#118) and 
E6 is not detected. The Table 3 offsets are estimated from 
the positions of isolated stars around IRS 13E published 
by Ott et al. (1999) which, in this paper are given relative 
to IRS 7. The offset of IRS 7 from SgrA* estimated from 
VLA measurements, taken from Menten et al. (1997) is 
added. The identification of the sources is presented on 
Fig. 2. 

Table 3. Offsets from SgrA*" of the sources shown on Fig. 2. 



ID 


ID'" 


Aacos5{") 


Asn 


W 


40 


-4.17 


-1.96 


El 


25 


-3.02 


-1.58 


E2 


27 


-3.24 


-1.68 


E4 


77 


-3.26 


-1.35 


5 


69 


-3.64 


-1.20 


6 


101 


-2.90 


-2.77 


7 


114 


-2.74 


-2.05 


8 


120 


-4.45 


-1.53 


E3A 


118 


-3.24 


-1.48 


10 


145 


-4.35 


-1.19 


E3B 


118 


-3.30 


-1.44 


12 


188 


-4.44 


-2.60 


13 


296 


-2.56 


-0.79 


14 


295 


-2.39 


-0.88 


15 


328 


-3.18 


-0.63 


16 


381 


-3.44 


-0.67 


E5 


780 


-3.46 


-1.50 


18 


902 


-2.55 


-1.42 


E6 




-3.10 


-1.37 


20 


459 


-4.15 


-0.85 


X 




-3.44 


-1.79 



a: SgrA* position (Yusef-Zadeh ct al. 1999): 
a(2000) = 17 45 40.0383 ± 0.0007s 
(5(2000) = - 29 00 28.069 ± 0.014" 
b: identification number from Ott et aJ. (2003) 




Fig. 1. IRS 13 field from the Gemini AO image in the Kp band. 
IRS 13E is the central, compact group of stars and IRS 13W the 
brightest source ~ 1" southwest of IRS 13E. The coordinates 
are in arcsec offset from SgrA*. 



Table 4. Projected separation of the main IRS 13E sources. 
From the astrometry the precision on the distance betwen two 
sources is ± 15 AU. 



Sources 


Ang. Sep. 


Distance (AU) 


El - E2 


0.241" 


1930 


El - E4 


0.330" 


2630 


E2 - E4 


0.327" 


2600 


E2 - E3A 


0.200" 


1600 


E3A - E3B 


0.073" 


590 



We used this astrometry to examine the position of 
the X-ray source CXOGC J174539. 7-290029 whose IRS 13 
is proposed as the optical counterpart by Baganoff et al. 
(2003). The X-ray source was placed on Fig. 2 by estimat- 
ing the offset of the source with respect to SgrA* from the 
coordinates of both sources reported in Baganoff et al. . 
The l-cr error reported on the coordinates of each X-ray 
source is of ±0.2" in right ascension and ±0.1" in dec- 
lination. By combining the astrometric uncertainties on 
the two source positions we drew a minimum error box 
of ±0.3" around the nominal position. The resulting po- 
sition falls outside the center of IRS 13E, which will be 
discussed in Sect. 6. 



Fig. 2. The Fig. 1 field with the star detected after deconvolu- 
tion by the MCS code (Sect. 2.3). The vector associated with 
most of the stars represents in amplitude the velocity and the 
direction of proper motions measured from SHARP data by 
Ott et al. (2003). For E3A and E3B, only the proper motion of 
the center of light is determined. The amplitudes reported in 
Table 6 for the four brightest sources of IRS 13E give the scale 
of the proper motion vectors. The cross marked X represents 
the nominal position of the X-ray source at the center of an 
error box of ±0.3" (see Sect. 3.2) . 




Fig. 3. At left, NICMOS images of IRS 13E in the Paa fil- 
ter (1.87/im), and at right, difference between this image and 
the suitably scaled image in the continuum at 1.90 /im. Only 
IRS 13E2 and IRS 13E4 remain which appear as emission line 
stars. 



3.3. Photometry 

Table 2 has been used to derive the standard photome- 
try in H, K, and L bands and the color indices of all the 
sources identified in Fig. 2. The Kp filter (Table 1) which 
is close and has a width comparable to K was used. It is 
referred to as K in the results presented in Table 5. 
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Table 5. H, K and L photometry of the IRS 13E cluster and 
of the nearby field stars from Table 2. 



ID 


H 


K 


L 


H-K 


K-L 


W 


14.55 


11.30 


8.92 


3.25 


2.38 


El 


12.74 


10.98 


8.59 


1.76 


2.39 


E2 


13.05 


11.03 


7.73 


2.02 


3.30 


E4 


14.37 


11.72 




2.65 




5 


14.32 


11.90 




2.42 




6 


14.71 


12.13 




2.58 




7 


14.41 


12.18 




2.23 




8 


14.86 


12.23 




2.63 




E3A 


16.48 


12.47 


7.50 


4.01 


4.97 


10 


15.55 


12.80 




2.75 




E3B 


16.95 


13.07 


7.92 


3.90 


5.14 


12 


15.92 


13.28 




2.65 




13 


15.95 


13.35 




2.60 




14 


15.95 


13.42 




2.53 




15 


16.65 


14.20 




2.45 




16 


16.66 


14.22 




2.44 




E5 


18.19 


14.28 


8.48 


3.90 


5.79 


18 


16.79 


14.52 




2.29 




E6 


16.24 


14.56 




1.68 




20 


17.29 


14.78 




2.52 





3.4. Proper motions 

From the results communicated by Ott et al. (2003), the 
direction and amplitude of the proper motions for the 
main IRS 13E sources and most of the stars in the 2.5" 
field are represented in Fig. 2. Note that the five sources, 
13E1, 13E2, 13E3A, 13E3B (proper motion is only given 
for the center of light of 13E3A and B) and 13E4 arc all 
moving West, with a similar velocity, while all the nearby 
stars have very different directions and amplitudes. The 
amplitude of the proper motions of the four main sources, 
in angular motion per year (for a GC distance of 8kpc) 
and in velocity, are presented in Table 6. The reported un- 
certainties come from the mean error on the proper motion 
vector coordinates adjusted on the set of positions. The 
rms uncertainty on 13E3 proper motion is the largest one, 
because of the difficulty of measuring accurate positions 
from several epochs of SHARP data for a weak source so 
close to much brighter sources, 13E2 and 13E4. The error 
becomes large for data recorded with poor seeing condi- 
tions. 

Table 6. Amplitude of proper motions of the main IRS 13E 
sources (from Ott et al. 2003). 



Name 


mas/yr 


km/s 


13E1 


5.50 


207 ±14 


13E2 


8.20 


310 ±19 


13E3 


7.54 


285 ±270 


13E4 


7.77 


294 ±32 



3.5. The ionized gas near IRS 13E 

The Paa emission line is tracing also in the two central 
parsecs the HII region called SgrA West or commonly 
the Minispiral. IRS 13E is located just at the northern 
end of a bright emission arc, part of the circular struc- 
ture open toward SgrA*, called the Minicavity. In Stolovy 
ct al. (1999) the bubble- like feature is proposed as cre- 
ated by the wind from a star identified by a weak Paa 
point source, since located very near the geometric cen- 
ter of the Minicavity. A proper motion of the edge of 
the Minicavity of ~ 200 km s~^ is measured by Zhao & 
Goss (1998). In Paumard et al. (2003a), it is shown that 
the Minicavity is embedded inside a non-planar gas flow 
of the Minispiral, called the Northern Arm. In the field 
of the Minicavity another velocity structure is identified, 
stretched toward northwest, called the Bar, roughly per- 
pendicular to the bright edge of the Northern Arm. On 
the line of sight of IRS 13E these two velocity structures 
produce the two main components in the line profile of 
the ionized gas shown on Fig. 4 central panel. The fastest 
component at -250 km is due to the motion of the 
gas disturbed by the Minicavity, and the slowest compo- 
nent, at -39 km s~^, to the Bar. On a morphological basis 
IRS 13E would seem to belong to the bright arc of the gas 
shocked at the edge of the Minicavity. On the intensity 
map of the Bar, isolated by this multi-component analy- 
sis, a small region just centered on the IRS 13E position 
appears locally enhanced. Therefore, the observed bright- 
ness of the ionized gas around the position of IRS 13E 
is due to the addition of two contributions on the line of 
sight: the edge of the Minicavity and the locally excited 
gas of the Bar. Only this latter component is due to the 
presence of the IRS 13E sources, by their strong ionizing 
flux. Hence, IRS13E should be located close to or inside 
the Bar, which lays behind the Minicavity (Paumard et 
al. 2003a). 

3.6. The stellar Paa emission line 

As shown on Fig. 3, only two of the IRS 13E compo- 
nents, IRS 13E2 and 13E4, remain after subtraction of 
the F190N continuum from the F187N filter which con- 
tains Paa. These two stars arc unambiguously emission 
line objects, with the integrated line intensity at IRS 13E2 
brighter by a factor 2.35 ± 0.1 than at IRS 13E4. 

3.7. The Brj and the Hei 2.058 fim line profiles 

The profiles of these two emission lines were obtained from 
spectra extracted on the same aperture size (3x3 pixels 
i.e. ~ 1" X 1") from each BEAR data cube, at the IRS 13E 
position. With the spatial resolution of the BEAR data 
the contribution of the two emission line stars 13E2 and 
13E4 is mixed, to which is added the ISM emission. After 
the most plausible subtraction to each line profile of this 
latter contribution, shown on the central panel of Fig. 4, 
possible thanks to the spectral resolution of the data, Br7 



Maillard et ah: GC source IRS 13 at high resolution 



8 

appears much narrower than the He i 2.06 /xm hne. With a 
FWHM of 215 km s-\ compared to ~ 900 km s"^ for the 
helium Unc, as ahcady measured in Paumard ct al. (2001), 
the two hues should belong to two different sources, which 
is discussed in Sect. 4.2. 



Hel 2.06/21X 1 ISM Br7 2.16um 




Fig. 4. He I 2.06/um and 617 emission line profiles at IRS 13E 
normalized to 1, from the BEAR data, after subtraction of the 

interstellar emission (ISM) contribution the observed profiles. 
The diflerence of linewidth is noticeable. The multi-component 
ISM emission, extracted from a ring around IRS 13E in the Br7 
data cube, is shown in the middle panel. Since Hei 2.06/im is 
coming only from 13E2 (Clenet et al. 2003) the broad profile 
adjusted on the wings of the Br7 profile, of same width as 
the He l profile, shows the relative intensity of Br7 from 13E2 
compared to 13E4. 



3.8. Spectral energy distribution and extinction law 

The photometric measurements presented in Table 2 were 
used to build the SED between 1 and 4 /im of the IRS 13 
components (Fig. 5), and also of all the stars detected 
in the surrounding field. To achieve the final goal of de- 
termining the spectral type of these stars, a dereddening 
has to be applied over this range. The extinction over the 
central parsecs is highly variable (Blum ct al. 1996, Rieke 
1999, Scoville et al. 2003). We took the most recently pub- 
lished law, which is a merging of previous works (Moneti 
et al. 2001), and adjusted Ay, making first the simplifying 
assumption that the reddening toward IRS 13E and the 
few arcsecs around would be identical for all the sources. 
For this adjustment we were helped by two constraints: 
IRS 13E2 and 13E4 are emission line sources (Fig. 3), 
therefore hot sources with a blue SED, while IRS 13W is a 
cool star (Krabbe et al. 1995). Black-body curves were fit- 
ted to the data. Ay was adjusted in order to fulfill the two 
constraints. Higher values of Ay make the sources bluer, 
while lower values make the sources redder. Temperatures 
higher than 25,000 K had to be introduced for IRS 13E1, 
13E2 and 13E4, which correspond to the Rayleigh- Jeans 
regime in this spectral range. In that case, only a lower 
limit of Teff can be derived. The slope of the SED be- 
comes constant in a log[F(A)] diagram. For these stars the 
adjustment of Ay makes it possible to bring the SED par- 



allel to the data points, providing the strongest constraint 
on Ay. Finally, a value of = 35 ± 0.5 was adopted from 
the three hot stars, 13E1, 13E2 and 13E4. A mean value 
of Ay = 30.5 had been determined by Rieke (1999) from 
a survey of the stars in the central parsec, which excluded 
objects like IRS 13 from the color-magnitude diagram be- 
cause of the surrounding dust. We confirm a much higher 
Ay value for IRS 13E. On the other hand, for IRS 13E6 
Ay = 35 was too high for a good fit to the data. A value 
of 29 was more appropriate. On the contrary, we will show 
in Sect. 4.7 that a value of Ay = 38.5 had to be adopted 
for IRS 13W. 

However, it appeared that the fit of the dereddened 
data from Table 2 was not possible with a single temper- 
ature for the IRS 13 sources. Except 13E6, they show an 
infrared excess. Then, the fit was made as a sum of two 
black-body curves, Coe/i x BB{Ti) + Coe/2 x BB{T2). 
The four parameters of adjustment obtained for each star 
are reported in Table 7. The final SEDs are shown on 
Fig. 5. From all the results presented above a spectral 
classification of the stars detected in the IRS 13E cluster 
is proposed. 



10 





1 

+ 


w 








E1 









E2 






A 


E4 






□ 


E3A 






X 


E3B 









E5 






• 


E6 




^\ + " 

+ ^"^"^-^^^"-^ 






n 






^jiT- 


— -f^ - 


\ 











1 






1.0 1.5 2.0 


2.5 


3.0 


3.5 


microns 







Fig. 5. Dereddened SED of the IRS 13 sources in 
Wcm~^ IJiVD.~^ . The top of the arrows indicates the estimated 
upper limit of the detectable flux in the 1.1 /um and the L-band 
filters. The various lines represent the best fit between 1 and 
4 fxrci of the data points from a two-component model with the 
parameters of Table 7 and the most adapted values. 

4. Nature of the IRS 13 sources 

4.1. IRS 13E1 

This source is characterized by: 
- no emission detected in Paa (Fig. 3), 
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Table 7. Fitting parameters of the SED of the IRS 13 sources. 
The spectral type of each source, as discussed in Sect. 4, is 
summarized in the last column. 



Star 


Coefi 


Ti K 


Cocf2 


T2 K 


Au 


Sp. Type 


W 


13.50 


3600 


4250 


640 


38.5 


M3III 


El 


0.800 > 


25000 


12000 


550 


35 


~ 051 


E2 


0.620 > 


25000 


40000 


550 


35 


- WN8 


E4 


0.140> 25000 


45 


1550 


35 


~ 05IIIe 


E3A 


0.460 


3800 


33000 


610 


35 


d. WR" 


E3B 


0.375 


3800 


29000 


580 


35 


d. WR" 


E5 


0.070 


6000 


9800 


630 


35 


d. WR" 


E6 


0.013> 


25000 






29 


~ 05V 



a : dusty Wolf-Rayet star 



- an elongated envelope in H and Kp toward 13E2, which 
shows at the star position a FWHM of 180 mas (to 
compare to the 40 mas of the PSF) or 1435 AU, 

- a SED fitted by a > 25,000 K source and an infrared 
excess at 550 K (Table 7), 

- a Kmag of 10.9 (Table 5). At the distance of the GC and 
with the adopted extinction, IRS 13E1 is too luminous 
to be a dwarf star. A 05V star of comparable Te// 
would have a Kmag — '^^: (with Ak^4:). IRS 13E1 is 
~ 3 magnitudes brighter. 

Hence, IRS 13E1 must be a blue, supergiant star. The 
surrounding halo of scattered emission is indicative of a 
strong stellar wind, which also favors the identification 
of a massive, hot star. Consequently, IRS 13E1 might be 
close to a 051 spectral type. However, at the position of 
the source there is no evidence in the Paa image of a neg- 
ative level (Fig. 3), which would be the signature of the 
photospheric hydrogen line in absorption since no emis- 
sion is detected. Possibly, with the uncertainty on the ad- 
justment of the continuum subtraction from the 1.87 /xm 
filter, the integrated emission from the envelope in the 
image can be compensating the photospheric absorption. 
Only a near-infrared spectrum of IRS 13E1 behind adap- 
tive optics could confirm this assumption. 

4.2. IRS 13E2 and 13E4 

IRS 13E2 and 13E4 are two emission line stars (Fig. 3). In 
the K band IRS 13E2 is brighter than IRS 13E4 by a factor 
2i 2 (Table 5). In our L band 13E4 is no longer detectable 
while 13E2 can still be measured. With this detection it 
turns out that all the emission spectra of IRS 13E obtained 
up to now (see references in Introduction) arc in fact the 
combination of two emission line stars which are likely 
different. 

A decisive element is reported by Clcnct ct al. (2003). 
From spectro-imaging with a FP in the 2.06 /xm helium 
line behind the CFHT-AO system they indicate that 
IRS 13E2 is the only helium emitter. The radical differ- 
ence of linewidth between the Hei 2.06 /im and Br7 stel- 
lar line profiles at IRS 13E (Fig. 4) brings another piece 



of information. A similar difference of linewidth of the 
He I 2.06 /um lines is observed among the hot stars of the 
central cluster (Paumard ct al. 2001, 2003b), leading to 
the identification of two classes of massive stars. Applying 
this criterion to 13E2 and 13E4, and taking into account 
the detection of Hci 2.06 /im only at IRS 13E2, it can be 
concluded that IRS 13E2 as a strong helium emitter with 
broad line is a late-type WR star, probably of WN type 
(Figer ct al. 1999) since no detection of C ill or C iv typi- 
cal of WC type is reported from earlier K-band spectra. A 
Teff > 25,000 K estimated for the source (Table 7) is con- 
sistent with this identification. IRS 13E4 with also a Teff 
> 25,000 K (Table 7), but no hehum emission, source of 
the narrower Br7 profile, is more likely a less evolved star 
than IRS 13E2. From its SED and its absolute brightness 
in K, fainter than IRS 13E1, IRS 13E4 can be reasonably 
proposed as a 05IIIe star or just reaching the LBV stage 
but in a high extinction phase to explain its weakness. 

However, there is an apparent contradiction between 
the brightness of IRS 13E2 in Paa in Fig. 3 and the Br7 
profile at IRS 13E from the BEAR data after correction 
for the ISM emission (Fig. 4). As Paa appears strong at 
IRS 13E2 in Fig. 3, the Br7 profile at IRS 13E should be 
dominated by the emission from this source and be as wide 
as the He I 2.06 /xm profile. That is not what is observed 
as shown on Fig. 4. To reconcile these two facts, it must 
be noticed that Paa is far from a perfect indicator with 
narrow-band imaging technique, to distinguish between 
hydrogen-rich and helium-rich emitters since the Paa line 
(1.8751 /im) is blended with a strong helium line, Hei (4- 
3) at 1.8697/im. There is another helium line within the 
bandpass of the continuum filter, at 1.9089 /xm, but which 
is weak and will contribute to subtract only a little of 
the helium emission. All these features are well seen in 
the CGS4 spectrum of the AF star around Paa, presented 
by Najarro ct al. (1994), which is another helium star 
belonging to the class of the broad-line stars (Paumard 
et al. 2001, 2003b). Hence, the intensity in the F187N- 
F190N it the star position cannot be considered as 

a fully reliable measurement of the true Paa emission in 
the stellar atmosphere. The bright spot at IRS 13E2 in 
Fig. 3 is likely due to the Hei 1.8697/im line and with 
some contribution of the Paa emission. That is consistent 
with the Br7 profile shown on Fig. 4 which can be decom- 
posed into a narrow line, whose origin must be IRS 13E4, 
and a fainter, broad component of same width as the Hei 
2.06 /xm line, which should be the contribution of IRS 13E2 
to the observed profile. The residual Br7 profile will show 
a P Cyg profile, typical of hot stars with an atmosphere 
in expansion, consistent with the spectral type attributed 
to IRS E4. 

4.3. IRS 13E3A and B 

IRS 13E3 resolves into a double source in the deconvolved 
AO images in H and Kp (Fig. 2). Their projected sepa- 
ration is equal to ~600 AU (Table 3.2). The two sources 
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are photometrically quite identical. They are extremely 
red objects as indicated by the measurements in Table 5, 
and from their SED (Fig. 5). They arc faint in the H 
band and prominent in the L band. We measure a K-L 
color index of ~ 5 mag. for the two components (Table 5). 
Several sources in the inner parsec, mainly located along 
the Northern Arm (IRS IW, 2, 3, 5, lOW, 21), are also 
very red objects, with K-L > 3, reported in Clenet et al. 
(2001). IRS IW and IRS 21 have featureless spectra in 
the K band (Blum et al. 1996). IRS 21 has been stud- 
ied in detail, from 2 to 25 /xm by Tanner et al. (2002). 
They have fitted its SED by a two-component model, the 
near-infrared scattered light from the central source peak- 
ing at ~3.8/um (760 K), and the mid-infrared re-emitted 
light from the dust shell at ^250 K. They conclude that 
IRS 21 is a dusty WR star, experiencing rapid mass loss as 
well as the other luminous Northern Arm sources (Tanner 
et al. 2003). The IRS 13E3 SEDs are also fitted by two 
components (Table 7), but with respectively 3800 K and 
600 K. IRS 13E3A and E3B are likely sources of the same 
type as IRS 21 and the other Northern Arm sources. The 
higher temperature of the infrared component can come 
from the additional heating of the dust shell by the very 
close, massive blue stars, IRS 13E1, 13E2 and 13E4. That 
is also consistent with IRS 13 being not a prominent source 
at 12.5 fxm on the images in this band (Tanner et al. 2002) 
compared to the other Northern Arm sources. 



4.4. IRS 13E5 

One of the sources revealed by deconvolution of the AO 
images (Fig. 2) we propose to name IRS 13E5, is also 
present in the SHARP data (Ott et al. 2003). From the 
dereddened photometry (Fig. 5) this source has a SED 
similar to IRS 13E3A and E3B, being roughly a factor 2 
fainter than each of the IRS 13E3 components. Its SED 
is also fitted by two thermal components, with tempera- 
tures of the same order as for 13E3A and E3B (Table 7). 
From this similarity we propose that IRS 13E5 is another 
example of dusty WR star, possibly more embedded, be- 
hind IRS 13E3A and E3B. 



4.5. IRS 13E6 

IRS 13E6 is detected in H, in the F160W filter, a broader 
H filter, and near the detection limit of the Kp band 
(Table 2) and not in L. It is not detected by Ott ct al. 
(2003) whose data come from K-band imaging on the 
NTT. With Teff > 25,000 K and K^ag = 14.56 (Table 5) 
13E6 is a weak, hot star, 3.5 magnitudes fainter than 
IRS 13E1. It can be considered close to a main sequence 
05V star. On Fig. 5 the best fit to the data was ob- 
tained with a value of Ay lower than for the other IRS 13E 
sources. It should mean that IRS 13E6 is located on the 
line of sight, but in front of the IRS 13E complex. 



4.6. The red halo around IRS 13E 

The flux in the L band is only due to continuum emission 
since the Bra line from the ionized gas falls outside the fil- 
ter bandpass. The deconvolution of the AO image in this 
band, after subtraction of the stars, leaves an enhanced 
background emission around IRS 13E shown on Fig. 6, 
extending on 2"in the North-South direction, above an al- 
most uniform continuum. On a larger scale, over the cen- 
tral parsecs, a diffuse emission is seen in the mid-infrared 
(Tanner ct al. 2002) whose structure is exactly following 
the intensity map of the ionized gas forming the Minispiral 
(Stolovy et al. 1999, Morris & Maillard 2000, Paumard et 
al. 2003a). This emission is the thermal emission from the 
dust dragged by the same fast flowing motion as the gas, 
heated through the trapping of Lya photons from the cen- 
tral, massive star cluster (Ricke ct al. 1989). In the analy- 
sis of the local structure of the interstellar gas (Sect. 3.5) 
we have shown that the IRS 13 complex is embedded in 
the flow called the Bar (Paumard ct al. 2003a). This flow 
is also likely mixed with dust. Hence, the enhanced contin- 
uum emission seen in L around IRS 13E should come from 
the dust heated up by the strong UV field coming from 
the local concentration of hot stars. There, the emission 
temperature of the dust has a value around 600 K from 
Table 7, against 250 K at IRS 21 (Tanner ct al. 2002) 
and probably lower out of the embedded sources. From 
further observations with NAOS/CONICA in the same 
band (Eckart ct al. 2003) the bright extension 1" north 
of IRS 13E resolves into a compact cluster of several very 
red sources. Their spectral type might be comparable to 
the dusty WRs as 13E3A, B and 13E5 detected next to it. 
Spectroscopic studies are needed to confirm this hypoth- 
esis. 

4.7. IRS 13W 

This star is associated with IRS 13 only for historical rea- 
sons. We kept it in all the study since its cool stellar type, 
confirmed by the presence of CO in its K-band spectrum 
(Krabbe ct al. 1995), was a constraint on the determina- 
tion of the local value of A-^ . In the deconvolved Kp-band 
image 13W shows a larger size than the PSF which might 
be indicative of a dusty envelope, attested also by a sig- 
nificant infrared excess (Table 7). However, in the first 
analysis with the assumption of the same A^ value for 
all the IRS 13 sources a Tg// of ~ 2600 K was obtained 
for 13W. With this temperature the cool star should be 
of Mira-type. A H-l-K spectrum would show the deep ab- 
sorption of water vapor on each side of the K band. With 
the NICMOS data we have a photometric measurement 
at 1.90 /xm, just in the water vapor band. Fig. 5 shows a 
smooth distribution of data points. Thus, we must con- 
clude that this temperature is too low. To increase this 
temperature in the fit a higher Ay value must be assumed 
locally. With the four parametc;rs of the fit (Table 7) plus 
Ay the solution is not unique. A plausible solution is ob- 
tained with Ay =38.5 and T^ff = 3600 K. A higher tem- 
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-2.0-2.5-3.0-3.5-4.0-4.5 
dacosS (arcsec) 

Fig. 6. Map of the residual fitted continuum in the deconvolu- 

tion operation of the L-band imago with the MCS code. The 
identified stars in this band (Table 2) have been subtracted 
explaining the dark holes in the map (e.g. at the W position) . 
Source positions are marked by crosses with their names. To 
limit confusion, for the E3 binary only E3B position is indi- 
cated. The apparent halo might be not purely thermal emission 
of dust and contain more fainter embedded sources (Eckart et 
al. 2003). An estimation of the halo contrast is given for the 
brightest pixel (position — 3.1", — 0.8"), which is four times 
brighter than east of 13E, out of the halo. 



peraturc would require a higher Ay vahie which would 
become inconsistent with the non detection of IRS 13W 
at 1.1 y^m (Table 2). With this Tg// temperature and its 
luminosity IRS 13W can be assumed to be a M3 giant 
star. 

4.8. The other stars of the IRS 13 field 

As presented in Table 2 twelve more stars are detected 
in the small field around IRS 13. Note that they are only 
detected in the 1 - 2.3 ^m range. They offer an opportu- 
nity to study in detail a small sample of the central parsec 
stellar population. By fitting black-body curves (a single 
temperature component was required), under the assump- 
tion of the same Ay as for IRS 13E, nine stars are main 
sequence K and M stars, with Tg// from 2800 to 5000 K. 
Hence, they belong to the numerous old star population of 
the GC (Blum et al. 1996). The three remaining stars (# 7, 
12, 18 - Fig. 2) are very hot stars, with T^ff > 25,000 K 
and Kmag respectively of 12.2, 13.3 and 14.5 (Table 5). 
As proposed for IRS 13E6, they should be main sequence 
O stars. Near infrared spectroscopy is needed to better 
constrain their spectral type. However, in this small field 
four out of a total of 14 sources (putting the massive com- 



ponents of IRS 13E apart) are apparently 0-type stars, 
which could make a proportion of about one quarter of 
such stars in the central parsec, the rest being late-type 
stars. 

5. IRS 13E as the remaining core of a massive 
star cluster 

IRS 13E appears as only composed of hot, massive 
stars. This concentration within 0.5" cannot be fortuitous. 
Further deep high-spatial resolution imaging could even- 
tually reveal more components (Eckart ct al. 2003). The 
common westward direction and similar amplitude of the 
proper motions with a mean value of ^ 280 km s~-^ for the 
main components (Sect. 3.4 and Table 6) is a decisive ar- 
gument to indicate that 13E1, 13E2, 13E3A/B and 13E4 
are physically bound. Regarding 13E5 for which a proper 
motion cannot be currently measured to associate it un- 
ambiguously to the cluster, its spectral type, identical to 
13E3A/B, leads us to conclude that this source belongs 
to IRS 13E too. For IRS 13E6 the proper motion is not 
available either, but the value of Ay 6 magnitudes lower 
than for IRS 13E likely indicates that the source does 
not belong to the cluster. The source is another 0-type 
star, just located on the line of sight. However, in conclu- 
sion, the source historically called IRS 13E is a compact, 
massive star cluster with at least six members. That also 
means a young star cluster of a few 10^ yr old, since sev- 
eral members are identified as having already reached the 
WR stage. 

The compactness of the cluster and the common 
proper motion of the components raise the question of the 
force which keeps the massive stars bound. The hypothesis 
of the presence of a dark, massive source, a stellar black 
hole at the center of IRS 13E is natural. An intermediate- 
mass black hole (IMBH 10^ to lO'' Mq) is supposed to 
form by runaway growth in massive, young stellar clus- 
ters as a result of stellar collisions in the cluster center, 
as was modeled by Portegies Zwart & McMillan (2002). 
Constraints on the possible central mass could be obtained 
from the radial velocities of the sources. It is estimated 
only for the two emission line stars, with a positive veloc- 
ity of ~ 30km s~^ for IRS 13E2 and a negative velocity of 
^ 30 km s-i for IRS 13E4 (Fig. 4). Associated with paral- 
lel and equal proper motion vectors (Fig. 2), this suggests 
that both stars orbit around the center of mass in a plane 
orthogonal to the plane of the sky. Assuming that both 
stars orbit around the black hole in a symmetrical fashion 
on circular orbits, half the projected separation between 
IRS 13E2 and IRS 13E4 (Table 3.2) gives an orbit radius 
of 1300 AU/ cosi where i is the angle that the line con- 
taining both stars makes with the plane of the sky. The 
radial velocity of each star (~ 30 km s~^) can be used as 
an estimate for their orbital velocity. A period of rotation 
of ~1295yr is derived. Then, the 3rd Kepler's law gives 
directly a total mass of ^ 1300 Mq cos~^{i)Ro/8 kpc. We 
can also release the circular orbits hypothesis, and only 
assume that the system is bound, which means that the 
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potential energy of each stars is greater than its kinetic en- 
ergy. This gives a lower limit to the black hole mass, half 
the previous estimate: Mbh > 750 Mq cos~'^(i)i?o/8 kpc. 
The dependence of these two values is quadratic on the 
orbital velocity of the stars and linear on their distance 
to the black hole, so that the constraint Mbh > lO'' Mq 
can be considered rather robust. It falls within the range 
derived by Portegies Zwart & McMillan (2002) for a black 
hole formed in the core of a dense star cluster with massive 
stars of initial mass > 50 Mq. 

5.1. The fate of a star cluster near SgrA* 

The large number of massive stars in the central parsec, 
which are very rare elsewhere in the Galaxy, remains one 
of the major mysteries of this region. Since star formation 
would be difficult due to the strong tidal forces from the 
SgrA* black hole, Gerhard (2001) made the interesting hy- 
pothesis that the central parsec He i stars, the most promi- 
nent of the massive young stars, might be the remains of 
a dissolved, young cluster, which originally formed fur- 
ther away from SgrA* . He argued that the Arches and the 
Quintuplet clusters, located within a projected distance of 
~ 30pc from SgrA*, testify that star formation by cluster 
of massive stars has been occurring in the nuclear disk of 
the Galaxy. If one examines Table 7, IRS 13E appears as 
a kind of summary of all the spectral types of young stars 
observed in the central parsec from O to WR. Among the 
helium stars only IRS 13E appeared multiple, which moti- 
vated the current study. With the hypothesis that IRS 13E 
might be the remaining core of a massive star cluster, was 
this cluster the source of the population of massive stars 
observed in the central parsec? 

Morris (1993) has argued that it would take longer 
than the lifetime of massive young stars to transport them 
inward within the central parsec if they formed at too large 
distance. The same argiimcnt is applied again by Figcr 
et al. (2000), who claim that the SgrA* cluster (Genzcl 
et al. 1996, Ghcz ct al. 1999, Gezari et al. 2002) could 
not have formed more than 0.1 pc from the center, and 
then, the initial clump should have an exceptional density 
>410^^cm'\ Gerhard (2001) discarded this argument by 
a revision of the conditions for a cluster formed at 30 pc 
to spiral into the center within the lifetime of its most 
massive stars. The main condition is that the initial mass 
of the parent cloud in which the cluster forms must be 
massive enough (~2xlO^M0) to survive the evaporation 
in the strong tidal field of the nuclear bulge. He concludes 
that clusters significantly more massive than the Arches 
cluster and formed a little closer than 30 pc can reach the 
central parsec in due time. In order to test this statement 
Kim & Morris (2003) have made several simulations, for 
different masses (IO^Mq and IO^Mq) and different initial 
orbit radii (2.5 to 30 pc), of the dynamical friction on a 
star cluster near SgrA*. They came to the conclusion that 
some simulations can be regarded as candidates for the 
origin of the central parsec cluster, but that "the required 



conditions are extreme" , with an initial mass of the clus- 
ter of IO^Mq or a very dense core > IO^MqPC"^ (Kim 
ct al. 2003). A mass of the cluster of lO^M© supposes a 
very large number of particles, with some of initial mass 
> IO^Mq to reach rapidly the WR stage. Compared to 
the relatively small number of detected helium stars, con- 
centrated in the central parsec, ~ 19 from the revision by 
Paumard et al. (2003b), plus few more dusty WRs, a very 
large quantity of 0-typc stars (> 10^) should be detected, 
which is not the case from the proportion of such stars we 
count in the IRS 13 test field (Sect. 4.8). This might seem 
an argument against this scenario. 

However, another analytic work by McMillan & 
Portegies Zwart (2003) has reconsidered the fate of a star 
cluster near the central dark mass. They tried to address 
the problem more completely by taking into account, in 
addition to the initial mass and the distance to the center, 
the original mass function of the cluster, the initial cluster 
radius and the stellar evolution through mass loss during 
the inspiral time of the cluster. They conclude that star 
clusters born with masses ~ lO^M© within 20 pc from the 
center, with a half-mass radii of ~ 0.2 pc can reach a final 
distance of 1 pc within lOMyr. As a secondary conclu- 
sion, they assess that from their mass and their distances, 
the Arches and the Quintuplet clusters, will never reach 
the vicinity of SgrA*. This latter work makes the origin 
of the central, massive star cluster by the dissolution of a 
compact cluster in the galactic tidal field more plausible, 
not requiring extreme mass conditions as in the simula- 
tion of Kim et al. (2003). In conclusion, we propose that 
the IRS 13E cluster, by its unique location and composi- 
tion, is the possible core of an earlier massive star cluster, 
formed about lOMyr ago, within 20 pc of SgrA*, with a 
mass of ~ 1O^M0, which was the progenitor of the entire 
hot star population, from WR to 0-type stars, observed 
today in the central parsecs of the Galaxy. In addition, the 
hypothesis of IRS 13E harboring a IMBH as consequences 
on this scenario. 

5.2. The need for a second black hole 

A recent paper by Hansen & Milosavljevic (2003) came 
to our attention when the current work on IRS 13 was 
completed. The authors try to solve the difliculty of bring- 
ing disparate groups of hot stars, the WR-type stars, the 
LBV-type stars and the S-cluster, at their observed lo- 
cation within the timescale required within a single evo- 
lutionary scenario. They argue that massive star clusters 
can sink within the required star lifetime but are tidally 
disrupted at a distance greater than 1 parsec from SgrA*, 
from which it would result a population of sources with low 
binding energy orbits unlike those of the helium stars and 
particularly of the S-cluster stars orbiting within 0.1 pc 
of SgrA**. For this purpose they propose a model with 
an infalling IMBH. The stellar orbits continue to evolve 
by undergoing close encounters with the IMBH, bringing 
some stars near enough to be trapped by the massive BH. 
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This model appears as a refinement of the same idea of 
an origin of the young star population in the central par- 
sec as formed at a distance of the SgrA* > 10 pc where 
star formation can occur, in a dense star cluster sinking 
towards SgrA*. We propose IRS 13E with a IMBH as the 
possible remnant of this initial cluster. 



6. What is the source of the X-ray emission at 
IRS 13? 

The origin of the X-ray emission at IRS 13 reported by 

Baganoff et al. (2003) for the position, by Coker ct al. 
(2002) for the spectrum, and by Muno et al. (2003) for 
complementary data on the energy distribution, must be 
examined in the context of the nature of IRS 13E pre- 
sented in the current work. According to Baganoff et al. 
the X-ray source is located 0.56" northeast of IRS 13W. 
Using our own astrometry (Table 3) and the offset of the 
IRS 13 X-ray source from the X-ray source at SgrA*, the 
source falls 0.75" east of IRS 13W (Fig. 2). With some 
minor difference we agree that the nominal position of the 
X-ray source does not coincide with the center of IRS 13E 
(our position is 0.37" south). This discrepancy is larger 
than the error box on the X-ray source. Then, the source is 
apparently not coincident with either IRS 13E2 as a post- 
LBV binary as originally proposed by Coker et al. (2002) , 
or the IRS 13E chister itself, which would have been the 
most plausible candidate with its rare concentration of 
several massive, hot stars with high mass loss and very fast 
winds, and possibly with a black hole at its center. Except 
a better astrometry the source CXOGC J174539.7-290029 
must be considered as a source independent of IRS 13E. A 
systematic X-ray survey of the Galactic Center region at 
sub-arcsecond scale with Chandra by Muno ct al. (2003) 
over a field of 17'xl7' centered on SgrA* has revealed 
more than 2000 discrete X-ray sources. Stellar remnants, 
white dwarfs with magnetically accreting disks, binaries 
with neutron stars or solar-mass black holes are consid- 
ered as responsible for a large fraction of these discrete 
sources. Such a source, too weak in the near infrared to 
be detected, could be present on the line of sight of IRS 13, 
since as measured by Muno et al. (2003) the distribution of 
the density of discrete sources peaks in the central parsec. 
The binary explanation for many stellar X-ray sources is 
based on a distinct emission feature centered at ~ 6.7 keV 
and otherwise a featureless spectrum. The spectrum pre- 
sented by Coker et al. (2002) stops at 6keV, probably 
because no energy was detected beyond, but it shows a 
narrow emission line at 3 keV. Wang et al. (2002) indicate 
that this line should be due to the contribution of mas- 
sive stars. Hence, the infrared counterpart to the X-ray 
source could be detected by deeper spectrophotometry of 
the IRS 13 field. A good S/N ratio X-ray spectrum of the 
so-called IRS 13 source would also help to better constrain 
the nature of this source. 



7. Conclusion 

The presence of a compact star cluster of six hot, mas- 
sive stars at the position of IRS 13E from high-resolution 
near-infrared observations is demonstrated. The spectral 
types of the various members range from O to WR, in- 
cluding dusty WRs. Proper motion measurements indi- 
cate that the brightest stars are co-moving suggesting that 
the members of the cluster are bound by a central IMBH 
with a mass >13OOM0. Such a secondary black- hole in 
the vicinity of SgrA* could be the element needed to ex- 
plain the population of massive young stars observed to- 
day (Hansen & Milosavljevic 2003). 

To precise the spectral type of the components, better 
constrain the mass of the IMBH, spectroscopy in the 1- 
5 /um range of all the individual sources within IRS 13E, 
at angular resolution as good as 0.1", is required. In the 
L and M bands it should confirm the expected feature- 
less spectrum, except dust signatures, of the IRS 13E3 
and 13E5 objects. These studies will need near-infrared 
3-D spectrometers behind an AO system on a 8-m tele- 
scope, like SINFONI (Mengel et al. 2000 and AMBER 
behind VLTI (Pctrov ct al. 2000). Deeper AO imaging, 
as already obtained with NAOS/CONICA north of the 
IRS 13E center (Eckart et al. 2003), could make it possi- 
ble to detect more members of the cluster. Proper motions 
of the fainter members would help to confirm which of 
the individual sources are kinematically bound together. 
Theoretical work is needed to confirm whether IRS 13E 
can be the remnant of a massive cluster. More generally, 
to address the problem of recent star formation in the 
vicinity of SgrA* a full census of the hot star popula- 
tion, illustrated in this paper on a small field, remains 
to be completed. It can be done by deep AO imaging to 
the condition to acquire data down to 1 //m. Up to now, 
such data, which appeared essential as shown in the cur- 
rent analysis, have only been possible with NICMOS on 
HST. Search for optical counterparts of the star-like X- 
ray sources detected by Chandra in the central parsec is 
another objective. 
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